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Higher order calculations in perturbative QCD have progressed in an extraordinary way in recent 
years, motivated by the needs of accurate phenomenology at the LHC. The frontier of present-day 
perturbative calculations is the next-to-next-to leading (NNLO) order [j|], due to the availability of a 
variety of novel computational techniques. Results for cross-sections at NNLO can be used thanks 
to the recent determination ^ of three-loop splitting functions which drive NNLO perturbative 
evolution. However, perturbative evolution at NNLO is unstable in the high energy (small x limit): 
the size of the NNLO corrections diverges as x — > at fixed scale. 

Small x resummation, which should take care of this instability, has a rather long history, 
starting with the original determination of leading high energy corrections ||3|] and their inclusion 
in perturbative anomalous dimensions [ffl]. Until quite recently, however, its relevance for phe- 
nomenology has been modest: since the advent of HERA data, it is clear that a NLO description 
of observed scaling violation is perfectly adequate [Q, ^], and the data show no evidence for large 
small x effects. The reason why nominally large corrections seem to have no impact has been 
obscure for a long time. However, due to some accidental zeros in coefficients, small x contribu- 
tions to NLO perturbative evolution are small, so in practice one could simply ignore the issue for 
all practical purposes. At NNLO, however, small x terms are large, the perturbative instability is 
manifest, and resummation becomes mandatory. 

Fortunately, over the last few years a fully resummed approach to perturbative evolution has 
been constructed. Within this approach, it is possible to understand why fixed perturbative order 
corrections are very large at small x, yet their full resummation leads to a considerable softening of 
small x terms, consistent with the fact that the HERA data do not show any large departure from 
NLO predictions. In order to obtain stable resummed results one must satisfy various physical con- 
straints, such as momentum conservation, renormalization group invariance and gluon exchange 
symmetry. These require the inclusion of several classes of terms which are formally subleading 
in comparison to the series of leading or next-to-leading small x logs. Once these constraints are 
enforced, the resummed perturbative expansion at small x becomes stable, and one no longer sees 
the strong small x enhancement or suppression that the leading [^J and subleading J7|] - JTo|] small x 
logs would give. 



A comparison of existing approaches to this resummation, as discussed respectively in refs. [ 1 1 



12|, (ABF approach) and [BJ] (CCSS approach), shows [@| that they yield results which 
agree with each other within the expected theoretical uncertainty. They agree in including the 
physical assumptions listed above. They differ mostly because the CCSS approach is built up 
within the BFKL framework, by improving the BFKL kernel through the inclusion of terms which 
become important in the collinear region, while the ABF approach is based on the construction of 
an improvement of the GLAP anomalous dimension through the inclusion of terms that become 
important in the small x region. The fact that they lead to very similar results is thus a consequence 
of the "duality" which relates the BFKL and GLAP description of perturbative evolution: at lead- 
ing twist, they can describe the same physics, provided the respective evolution kernels are suitably 



matched [O, ll7L [T^] 1 . 



Because the leading high energy corrections are dominated by gluon exchange, the resumma- 



1 Small-x corrections to polarized structure functions, of different form and origin than those considered here, have 
been recently discussed in ref. Jl6| ] 
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Figure 1: The splitting functions xP gg and xP qg for nt — 4 and a s = 0.2 as a function of x. Fixed order 
perturbation theory LO (black, dashed); NLO (black, solid) and NNLO (green, solid) and resummed LO 
(red, dashed) and NLO in QoMS scheme (red, solid) and in the MS scheme (blue, solid). 



tion is most easily performed in the pure Yang-Mills theory, and indeed the full resummed results 
for perturbative evolution of refs. [ pi] , p3] ] have been obtained with n f = 0. In order to actually 
get predictions for (flavour singlet) physical observables, one needs the full two-by-two matrix of 
splitting functions, and the set of hard coefficient functions for the desired observables, in a given 
factorization scheme. The construction of resummed splitting functions is simplified by the fact 
that only one of the two eigenvalues of the singlet anomalous dimension matrix has leading N = 
singularities: hence, only this eigenvalue is affected by the resummation. In deep-inelastic scatter- 
ing, the construction of resummed coefficient functions is further simplified by the fact that virtual 
photons at leading order only couple to quarks. This implies that there always exist schemes where 
only one partem (quark or gluon) contributes to each of the structure functions F2 and F^. It follows 
that in any factorization scheme the resummation of the coefficient is specified in terms of a single 
function for each structure function (these functions have been determined in ref. [19], where they 
are called hi and hi for F2 and Fl respectively.) 

Hence, (at least) two strategies are available for the construction of resummed observables in 
deep-inelastic scattering. The first possibility is to simply pick a factorization scheme, then con- 
struct resummed two-by-two evolution kernels and resummed coefficient functions in that scheme. 
This program was started in ref. p0[], where the full rif 7^ resummed evolution matrix was con- 
structed by extending a BFKL-like approach to coupled quark and gluon evolution, along the lines 
of the approach of refs. [|l4|, 15]. This has the advantage of giving evolution equations for off-shell, 
unintegrated parton distributions, but it has the shortcoming of providing results in a factorization 
scheme which only coincides with MS up to the next-to-leading fixed order, and differs from it 
at the resummed level. Available resummed coefficient functions [19], which are given in MS or 
DIS, are not readily combined with the evolution kernels determined in this way. Similar problems 
were encountered in ref. [21], where resummed structure functions were obtained by combining 
resummed anomalous dimensions and coefficients determined in different factorizaton schemes. 

A second possibility consists of taking advantage of the peculiar fact that both the resumma- 
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Figure 2: The /^-factors singlet Fi and Fi structure functions at fixed x = 10~ 2 , or as a function of 
Q with a s running and tif varied (the breaks in the curves correspond to the b and t quark thresholds). Fixed 
order perturbation theory NNLO (green, dashed); resummed NLO in QoMS scheme (red, solid), resummed 
NLO in the MS scheme (blue, dot-dashed). 



tion of the anomalous dimensions and that of the deep-inelastic coefficients for a given structure 
function are each determined in terms of a single function: it is only when specifying the factoriza- 
tion scheme that a matrix of anomalous dimensions and a vector of coefficients for each structure 
function are obtained. It follows that if one has full control of the transformation between factor- 
ization schemes at the resummed level, one can obtain results for resummed coefficient functions 
and anomalous dimensions in any scheme. 



This approach was pursued in ref. [22], where by making use of the general techniques for 
cross-section resummation developed in ref. [23] we were able to compute the full rif / matrix 
of resummed singlet splitting functions (see fig. 1) together with the coefficients C' q and C' g (i=2,L) 
for the singlet structure functions F2 and Fl, respectively, in various commonly used schemes such 
as the MS scheme. This then allowed us to determine for the first time the structure functions 



at a fully resummed level in a consistent factorization scheme. While referring to ref. [ |22| ] for a 
full discussion of the resummed splitting functions and coefficient functions, we focus here on the 
impact of the resummation on physical observables, and on the extraction of parton distributions 
from them. 

Thus, in fig. 2 we compare the effects of resummed evolution with those of fixed order pertur- 
bative evolution at NLO (which is the usual baseline for most current fits) and also at NNLO. We 
show the ^-factor (defined as resummed/(NLO fixed order pert)) for the singlet F2 structure func- 
tion at fixed values of x as a function of Q. For each x value we present three curves: the resummed 
case in the QoMS scheme, the corresponding plot in the MS scheme, and the NNLO fixed order 
perturbative. The starting point at <2o = 2 is taken to be the same for the structure functions plotted 
in the three curves. In the expanded linear scale of these plots we can see the moderate scheme 
difference between the QoMS and the MS schemes, which is left after combining coefficients and 
parton densities, which, individually, show a much more pronounced scheme dependence. 

For F2 the effect of resummation, at sufficiently small x values, goes in the opposite direction 
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Figure 3: The ^-factors, defined as the ratio of the fixed order NNLO or resummed to the NLO fixed order 
result, for xq and xg when F2 and Fi are fixed at the reference scale Qq — 5 GeV. Results are shown at 
Q = Qo as as function of x in the range x = 10~ 2 or 10~ 6 . The curves (top to bottom at small x) are: fixed 
order perturbation theory NNLO (green); resummed NLO in QoMS scheme (red), resummed NLO in the 
MS scheme (blue). 



to the NNLO perturbative evolution: the resummed i^-factor is less than one, corresponding to a 
smaller structure function at higher scales than with fixed order perturbative NLO evolution. The 
effect of the resummation is somewhat larger than that of the NNLO. This shows that for x ^ 0.2 
the inclusion of unresummed NNLO terms is actually counterproductive. 

In view of using parton distributions extracted from HERA data for physics at the LHC, it is 
also interesting to study how the quark and gluon distributions change when the resummation is 
switched on, while imposing that the measurable structure functions be unchanged. We do this 
in fig. 3, where we display the ^-factors for the input quark and gluon distributions as a function 
of x, with the input scale Qq=5 GeV. Here the structure functions F 2 and F L at the input scale are 
kept fixed, but the coefficient functions are changed depending on the perturbative approximation 
employed. Whereas at NNLO the quark and gluon distributions are enhanced with respect to the 
NLO, at the resummed level they are suppressed; the suppression becomes increasingly significant 
at smaller x. The effect on the parton distribution is sizable, but when evolving up the scale depen- 
dence tends to reduce this effect leading to the more moderate corrections displayed in figure 2. 
The general conclusion, however, is that if resummation effects are disregarded, the associated er- 
ror in extracting parton distributions at HERA and evolving them up at LHC is of order of about 
5% at x rsj 10~ 3 , and as large as 20% for low values of x ~ 10 -6 . 
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